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Abstract 

The MiniBooNE Collaboration reports a search for z/„ and disappearance in the Am 2 region 
of a few eV 2 . These measurements are important for constraining models with extra types of 
neutrinos, extra dimensions and CPT violation. Fits to the shape of the Uu and energy spectra 
reveal no evidence for disappearance at 90% confidence level (CL) in either mode. This is the first 
test of disappearance between Am 2 = 0.1 — 10 eV 2 . 
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Neutrino oscillations have been observed and confirmed at mass splittings (Am 2 ) of 
~ 10~ 5 eV 2 and ~ 10 -3 eV 2 , called the "solar" and "atmospheric" oscillations respec- 
tively. The observed mixing is consistent with three generations of neutrinos and a unitary 
mixing matrix. Complicating this picture, the LSND experiment observed an excess of V e 
in a beam indicating a possible third Am 2 around 1 eV 2 thus requiring more than 
three neutrino generations or other exotic physics. Recently, the MiniBooNE experiment 
[2I excluded two neutrino appearance-only oscillations (98% CL) as an explanation of the 
LSND excess if oscillations of neutrinos and antineutrinos are the same. 

Exotic physics models , including sterile neutrinos, extra dimensions, and 

CPT violation have been proposed to explain the LSND observation. Some of these models 
can also accommodate the MiniBooNE v e appearance oscillation results. These models are 
testable with measurements of v ^ and disappearance which constrain any non-standard 
oscillations of v^-^Vx- As described in this Letter, the MiniBooNE experiment has performed 
searches for and disappearance which probe a region of interest, Am 2 = 0.5 — 40 eV 2 , 
not covered by two previous disappearance experiments, CCFR [v^ and F M ) 7(] and CDHS 
(z/ M only) [8J]. Unless otherwise stated, all statements about neutrinos hold true also for 
antineutrinos. 

For the MiniBooNE experimental setup, the detector is located at a fixed distance from 
the neutrino source. In this case, disappearance due to oscillations has a distinct signature 
as a function of neutrino energy, because neutrinos with different energies oscillate with 
different probabilities for the same distance traveled. Disappearance would be observable 
either via a deficit of events (normalization) or, alternatively, via a distortion of the neutrino 
energy spectrum (shape), or both (normalization + shape). The absolute normalization 
uncertainties in a single detector experiment such as MiniBooNE are large, hence a shape- 
only disappearance fit is performed. The flux to the MiniBooNE detector is provided by 
the Fermilab Booster Neutrino Beam (BNB) which is produced by 8 GeV protons incident 
on a 1 cm diameter, 71 cm long (1.7 interaction length) beryllium target surrounded by 
a magnetic horn pulsed at 174 kA. The horn uses positive current to focus ir + and K + 
mesons for the neutrino mode sample and negative current to focus n~ and K~ for the 
antineutrino mode sample. The mesons that pass through a 60 cm diameter collimator 259 
cm downstream of the target decay in a 50 m long tunnel to produce the f M beam. The BNB 



flux {9} is determined using a GEANT4|l0j based beam simulation which has been further 
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modified to include updated p-Be particle production data. 



ll[ is 541 m. 



The distance from the proton interaction target to the MiniBooNE detector 
The MiniBooNE detector is a 12 m diameter spherical tank filled with 800 tons of mineral 
oil. The detector is separated into an inner region filled with 1280 inward facing 8 inch 
photomultiplier tubes (PMTs) and an optically isolated outer region used to reject cosmic- 
ray induced events. Charged particles produced in neutrino interactions emit primarily 
Cherenkov light, though a small amount of scintillation light is also produced. Light and 
^article production and propagation in the MiniBooNE detector is modeled using a GEANT3 



12j based simulation, which was tuned using MiniBooNE and external data. 



Neutrino interactions are simulated with the v3 NUANCE event generator [13J. Prior to 
selection, approximately 42% of all events in MiniBooNE are charged current quasi-elastic 
(CCQE) scattering and 22% are charged current single charged pion production (CCl7r + /~) 
in both neutrino and antineutrino mode. 

The search for oscillations is conducted with a sample of CCQE events because of the 
high statistics and purity. The reconstructed neutrino energy (E® E ) is calculated assuming 
the target nucleon is at rest: 

e qe = 2(M n - E B )E ll - {E 2 B - 2M n E B + AM + M\) 
v ' 2[(M n -E B )-E^+p li cos9 lx ] [ ' 

where AM = M% — M^, M indicates the muon, proton or neutron mass with appropriate 

subscripts, E B is the nucleon binding energy, E^p^) is the reconstructed muon energy 

(momentum) and 6* M is the reconstructed muon scattering angle with respect to the neutrino 

beam direction. A small correction is applied in both data and simulation to account for the 

biasing effects of Fermi-smearing. At 300 MeV, the muon energy resolution is 7% and the 

angular resolution is 5 degrees. The average E® E resolution is 11% for CCQE events \\\ - 

A CCQE event sample is selected by identifying a single muon in the detector and its 

associated decay electron using the same criteria as in the previous measurement of CCQE 

model parameters on carbon [14 1. Timing information from the PMTs allows the light 

produced by the initial neutrino interaction (first "subevent" ) to be separated from the light 

produced by the decay electron (second "subevent"). The timing and charge response of the 

PMTs is then used to reconstruct the position, kinetic energy and direction vector of the 

primary particle within each subevent. Exactly two subevents are required in the analysis 
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(the muon and its decay electron). To reject cosmic ray interactions, both subevents are 
required to have fewer than 6 veto-PMT hits. The first subevent must be in coincidence with 
a beam pulse, have a reconstructed track center less than 500cm, and greater than 200 inner 
tank PMT hits to eliminate electrons from stopped cosmic ray muon decays. The second 
subevent must have less than 200 inner PMT hits to be below the decay electron energy 
endpoint. Finally, the distance between the electron vertex and the muon track endpoint 
must be less than 100 cm, ensuring that the electron is associated with the muon track. 
This selection also applies to the antineutrino mode sample as the final state nucleon is not 
reconstructed and the detector does not distinguish muon charge. 

The selection yields 190,454 data events with < E® E < 1.9 GeV for 5.58 x 10 20 protons 
on target (POT) in the neutrino mode sample; 27,053 data events for 3.39 x 10 20 POT 
in the antineutrino mode sample. According to the simulation, the neutrino mode sample 
is 74% pure CCQE, and the antineutrino mode sample is 70% pure CCQE. The primary 
background (~75%) for both the and samples is CC17T events where the outgoing pion 
is unobserved (e.g. due to absorption in the nucleus). Though the neutrino mode sample 
has < 1% content, in antineutrino mode, the beam contains a substantial contribution 



of due to the higher 7r + production at the target and the higher v ^ cross section. The 
antineutrino mode is predicted to have 25% content. 

The CCQE cross section depends on the axial vector form factor, which is commonly 
assumed to have a dipole form as a function of four-momentum transfer (Q 2 ) with one 
adjustable parameter, Ma, the axial mass. Global fits to the world's neutrino scatte ring 
data on deuterium yield Ma = 1.015 GeV 15|, however recent results from K2K 16|, [17] 
and MiniBooNE 14J suggest a higher effective value of Ma ~ 1.2 for nuclear targets. In 
addition, MiniBooNE has also adjusted the level of Pauli blocking in the prediction, using a 
parameter k = 1.019, to better reproduce the experimental data at low Q 2 14J. The effect 
of Ma and k on the Q 2 shape is pronounced while oscillations provide relatively little Q 2 
distortion. Therefore, the effect of the cross section tuning does not mask any underlying 
disappearance in the neutrino or antineutrino mode samples. 

For the disappearance search, systematic uncertainties are included for the underlying 
neutrino flux prediction, neutrino interaction cross section, and detector response. The 
method used to estimate the uncertainties due to the underlying neutrino flux prediction and 
detector model is identical to that used in previous MiniBooNE results 0, [iij • The dominant 
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uncertainty on the CCQE cross section is from uncertainties on Ma and n, which span 
the difference between the deuterium and nuclear target results (Ma = 1.015 ± 0.20 GeV, 
k = 1.000 ± 0.019). In addition, the uncertainty on the shape of the CC1-7T background is 
estimated using the MiniBooNE CCl7r + data sample. With Ma = 1.015 GeV, the ratio of 
detected events to predicted events in MiniBooNE for neutrinos is 1.31 ± 0.26, 1.18 ± 0.18 
for antineutrinos, which shows agreement within the uncertainties. The difference between 
this value and previously published values is due to the different values of Ma and k 
Systematic uncertainties produce correlated errors between E® E bins that are included by 
developing, a covariance matrix in the same manner as in previous MiniBooNE oscillation 



analyses [2j, [18] . This covariance matrix includes separate normalization and shape-only 
error contributions. For the shape-only disappearance search, the prediction is normalized 
to data, and just the shape-only covariance matrix is used. 

The disappearance search uses the Pearson's x 2 test to determine allowed regions in the 
Am 2 — sin 2 29 plane. The x 2 is calculated from a comparison of the data, di, in E® E bin 
i to a prediction pi(Am 2 , sin 2 29) for 16 bins. The prediction assumes two-flavor — » u x 
disappearance characterized by one large mass splitting (Am 2 = Am^) between the light 
neutrino mass states k, which participate in standard three neutrino oscillation, and h, the 
heavier neutrino state, and one oscillation amplitude sin 2 29 = 4|E/ M) h| 2 (l — \U^^\ 2 ), where 
{U^hl is the muon flavor content of the heavy state h. 

16bins 

X 2 = E ( d i ~ NpjMy-^dj - Np 3 ) (2) 

where M^ is the shape-only error matrix, and iV is a factor which normalizes the prediction 
to the total number of observed events in data. All neutrino events in the prediction, 
including the CCl7r + background events, are allowed to oscillate in the fit based on the 
incident neutrino energy and distance traveled. The 90%CL limit corresponds to \ 2 > 23.5 
for 16 degrees of freedom (DF). The sensitivity is a fit to an unoscillated prediction including 
all statistical and systematic uncertainties. 

The top plot of Fig. [T] shows E® E after selection cuts for the neutrino data and the 
prediction assuming no oscillation (null hypothesis) with diagonal elements of the error 
matrix. The dominant systematics arise from the neutrino flux (production of 7r + /~ from 
p-Be interactions) and CCQE cross section; uncertainties at low energy are larger because 
of the substantial CCl7r + background and uncertainties in the CCQE cross section in this 
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FIG. 1: The top plot shows the E® E distribution for neutrino data (black) with statistical error 
rectangles (thickness of line indicates size of statistical error) , and prediction assuming no oscillation 
(grey). Attached to the prediction are the diagonal elements of the shape error matrix. The 
predicted CC17T background (dash) and background antineutrino (solid) events are also shown. 
The bottom plot shows the ratio of data to no oscillation (black), and the ratio of no oscillation 
to: Am 2 = 0.5 eV 2 , sin 2 26 = 1.0 disappearance (dashed line), Am 2 = 3.0 eV 2 , sin 2 26 = 0.5 
disappearance (dotted line) and for the minimum x = 

12.72 (13 DF) at Am 2 = 17.5 eV 2 , 

sin 2 26 = 0.16 (solid line). 

region. Though the diagonal elements of the error matrix are substantial, the correlations 
between energy bins are large. The x 2 between the data and the null hypothesis is 17.78 (16 
DF, 34% probability) for the neutrino mode sample which is consistent with no oscillation at 
90%CL. The top plot of Fig. [3]shows the 90% CL sensitivity and limit curves for the neutrino 
mode sample. The minimum x 2 = 12.72 (13 DF, 47% probability) at Am 2 = 17.5 eV 2 , 
sin 2 28 = 0.16, where the number of degrees of freedom is estimated from frequentist studies. 

The bottom plot in Fig. [T] shows the ratio of data to the null hypothesis and three 
oscillation scenarios. The shape distortion for Am 2 = 0.5 eV 2 is very different from Am 2 = 
3.0 eV 2 . The x 2 therefore changes rapidly as a function of Am 2 , resulting in rapid changes 
in the 90%CL sensitivity and limit curves (Fig. [3]) for small differences in Am? Similar 
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FIG. 2: Same convention as Fig. [T] for the antineutrino mode sample. Minimum \ 2 = 5.43 (11 
DF) is at Am 2 = 31.3 eV 2 , sin 2 29 = 0.96. 

features are also seen in previous disappearance analyses jz, 8]. 

The Pfj, disappearance analysis proceeds in the same manner as the analysis, except 
that only the events are allowed to oscillate in the fit and the events are kept fixed. 
This determines the limit on a model where the can oscillate but the cannot. A 
model where both and oscillate with equal oscillation probability versus energy would 
produce a limit very similar to the neutrino mode limit. 

During antineutrino data taking, two absorber plates inadvertently fell vertically into the 
decay volume at 25m and were later removed, creating three distinct data taking periods with 
zero, one, or two absorbers in the beamline. The event rate was predicted to be 13% (20%) 
lower for one (two) plate (s) in the beam. Approximately 15% of the antineutrino data taken 
had one absorber plate inserted, and 15% had two absorber plates inserted. Because the 
changes to the beamline are understood, a separate simulation was run with the appropriate 
number of absorber plates in the beamline. Figure [2] shows the E® E distribution for the 
antineutrino mode sample. The x 2 of the null hypothesis is 13.7, 8.2, 15.2, 10.29 (16 DF) for 
the zero, one, and two absorber plate and total data respectively. The antineutrino mode 
data is also consistent with no oscillation at 90%CL, so the bottom plot of Fig. [3] shows the 
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sin 2 (29) 

FIG. 3: The top plot shows the sensitivity (dashed line) and limit (solid line) for 90% CL for 
neutrino disappearance in MiniBooNE. Previous limits by CCFR (dark grey) and CDHS (light 
grey) are also shown. The bottom plot uses the same convention for antineutrino disappearance. 



90% CL sensitivity and limit curves for the antineutrino disappearance fit to all antineutrino 
data; the limit curves for the individual absorber data periods were found to be consistent 
with the total. In addition to the two-neutrino oscillation fits described above, some studies 
were performed some of the MiniBooNE energy spectra within the context of 3+2 oscillation 
models. The best fits for 3+2 sterile neutrino models in Ref. {4] are consistent with the 
MiniBooNE and data, but the data rules out the best fit point from the global fit 
to MiniBooNE u e data in Ref. at 90% CL with X 2 = 24.7(16 DF). 

In summary, MiniBooNE observes no evidence for or disappearance at 90%CL in the 
Am 2 region of a few eV 2 . The test of disappearance probes a region of Am 2 = 0.1 — 10 eV 2 
unexplored by previous experiments. 
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